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The Lag model has shown great promise in prediction of low speed and transonic sep- 
arations. The predictions of the model, along with other models (Spalart-Allmaras and 
Menter SST) are assessed for various high speed flowfields. In addition to skin friction and 
separation predictions, the prediction of heat transfer are compared among these models, 
and some fundamental building block flowfields, are investigated. 


I. Introduction 

One difficulty with current one and two-equation turbulence models is the inability to account directly 
for non-equilibrium effects such as those encountered in large pressure gradients involving separation and 
shockwaves. Current turbulence models such as Spalart’s one-equation model, 1 the classic k— e and Wilcox’s 
k — or two-equation models have been designed and tuned to accurately predict equilibrium flows such as 
zero-pressure gradient boundary- layers and free shear layers. Application in more complex flows can be 
problematical at best. Although there have been many attempts to modify or correct basic one- and two- 
equation models, most of these attempts have been only marginally successful in predicting complex flows. 

More complex models such as Reynolds stress models have been investigated extensively, primarily for 
relatively simple flows but also for complex flows. In most cases these models give somewhat better predic- 
tions than the simpler one and two equation models, but for complex flows they do not perform much better 
than the simpler models. One theoretical advantage of Reynolds stress models is that they directly account 
for non-equilibrium effects in the sense that the Reynolds stresses do not respond instantaneously to changes 
to the strain rate but more realistically lag them in time and/or space. Unfortunately, The Reynolds stress 
models are usually considerably more complicated and numerically stiff than the one- and two- equation 
models, and this has prevented their wide application for complex flows. 

In this paper we further investigate the Lag model, which was designed to account for non-equilibrium 
effects without invoking the full formalism of the Reynolds stress models. The basic idea is to take a baseline 
two-equation model and to couple it with a third (lag) equation to model the non-equilibrium effects for the 
eddy viscosity. The third equation is designed to predict the equilibrium eddy viscosity in equilibrium flows. 
One advantage of this method over comparable turbulence models is that it does not require wall distance, 
a major advantage for some flow solvers, and in complex flowfields where wall distance determination is a 
non-trivial exercise. 

Applications to four flows will be given including a high speed flat plate, cylinder flares, a rocket nozzle 
with significant separation, and space vehicles. 
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II. Method 


A. Reynolds Averaged Navier Stokes Equations 

The Reynolds averaged Navier-Stokes equations, written in conservation law form are 
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Turbulent kinetic energy, k, is simply assumed to be zero for the 1 equation turbulence models. 


B. Lag Model, Revised 

The Lag model has been revised from the original description. 1 It was simplified from the original description 
by dropping the leading function of Ry in the lag equation, and the constant defining the diffusion of 
the underlying k equation(ak) was increased to 1.5. This change was cosmetic, in that it rounded the 
turbulent/non-turbulent edge of the boundary layers and simplified up the model(the leading function of 
Rx was found to be unnecessary). All cases reported* earlier were re-run, and the previous results were 
reproduced, with minor improvement in the Johnson-Bachalo bump case. 

The revised Lag model is 
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with parameters 


do = 0.35 
a = 5/9 
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(3* = 0.09 


cik = 1.5 
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C. Numerical Method 


The code used in this study was 0VERFL0W2, modified to include the Lag model along with the high 
speed modifications involving the k discussed with the Navier Stokes Equations. Matrix dissipation was 
used with smoothing parameters as recomended by earlier studies of high speed flows with this code.' 1 
These values proved suitable for these flows also(Figl(a),l(b)), in that doubling and halving both second 
and fourth order smoothing coefficients did not affect the flowfields. The recommended eigenvalue limits 
(V £n = 0.3, V ei = 0.3) are not adjusted. 



(a) Boundary Layer 


(b) Skin Friction, x/L = .1 


Figure 1. Insensitivity of Solution to Smoothing Parameters, Flat Plate Solution 

The relaxation method is the implicit Pulliam Chausee diagonal method, with variable time step- 
ping(ITIME=l) or constant CFL (itime=3). CFL values for these high speed flows were chosen at 0.4, 
and the variable time step was adjusted as described in the overflow documentation." 

III. Results 


A. Flat Plate 

A Mach 5 flat plate" was simulated with baseline on a 129 x 129 grid(Fig. 2, with a Reynolds number based 
on (total)length of 10 x 10 s . The pressure wave created by the nose of the flat plate fits nicely within the 
grid, and wall spacing for this grid was constant at Ay/L = 2 x 10 -6 , which yeilded a Ay + ~ 0.5. 

The inflow and external face boundary conditions were freestream (plug flow) at M=5 and Too = 273.15, 
and the outflow boundary condition was simple extrapolation of all variables. Three cases were computed, 
one with an adiabatic wall by which the recovery factor of the model could be ascertained, and two fixed 
wall temperature cases, T w — 5.4Too and T w = 2.7 T^, one at roughly the adiabatic wall temperature and 
half that value. 

The solution was also obtained on two coarser grids, a medium 65 x 65 grid obtained by removing every 
other point in both directions, and a coarse 33 x 33 grid again removing every other point, this time from 
the 65 x 65 grid. These grids have Ay + of approximately unity and two, respectively. The grid stretching 
(Axj + i/Axj) was less than 1.06 in the streamwise direction, and below 1.08 in the wall normal direction, for 
the finest grid. The medium (65 x 65) grid had stretching below 1.12 streamwise and 1.16 wall normal. The 
coarsest grid had stretching below 1.24 streamwise and 1.35 wall normal. 

The Lag models skin friction predictions(here for the adiabatic wall) are in line with theoretical predic- 
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tions(Van Driest II with Karman-Schoener), and with other turbulence models(Fig. ??). The model shows 
similar wall spacing dependence as for the subsonic flat plate, with the skin friction predictions (Fig. 2(c)) 
requiring a Ay + < 1 for accurate skin friction determination. The boundary layer predictions (here plotted 
in wall coordinates, Fig 2(b)), are insensitive to this variation in grid density. Another comparison to an 
essentially flat plate flow will be seen in the inflow profiles and upstream heat transfer reported below for 
the cylinder-flare flowfield. 

Two skin friction predictions are shown for two wall temperatures, one essentially an adiabatic wall(Fig 
3(a)), and the other a cold wall(Fig. 3(b)). Skin friction is well predicted in both cases. The heating rate 
for the cold -wall is also well predicted(Fig 4(b)). The T w predicted by the various turbulence models for 
adiabatic conditions varies slightly(Fig 4(a)), with the recovery factor of the SA 1 model at 0.89 and Lag, 


SST and k— or model at .895. 



(a) Cf comparison, T w = 5.4T t T QW 



(b) Cf comparison, T w = 2.7Tt ~ T aw /2 


Figure 3. =5, Re^ =10x 10 s Flat Plate Skin Friction 


B. Cylinder Flare 

The Moo = 7 axisymmetric cylinder' flowfield tests the ability of the models to predict a simple well defined 
flowfield with a shock induced separation. The experimental geometries are an ogive-cylinder with cones of 
20°,30°,32.5o and 35° attached 1.39m from the nose of the ogive-cylinder. The cylinder-cone intersection is 
chosen as the origin of the x axis. The ogive nose is described more fully in, 5 ' and consists of a 10° conical nose, 
with a circular blend to the constant diameter cylinder body 0.644m downstream of the nose(x = — 0.746m). 
The ogive nose was modelled in an effort to match the experimental conditions as closely as possible. 

The finest grid for this case was 1179 x 257(streamwisexwall normal), extending from the ogive nose 
to the cylinder flare(Fig. 5). The interaction region 4mm ahead and behind the flare corner accounted for 
roughly 100 of the 1179 streamwise points. The portion upstream of this interaction region accounted for 
over 600 of the streamwise points, with the remaining points along the flare surface. Two sets of coarser grids 
were again created by removing every other point in both directions, resulting in a 589 x 129 and 295 x 65 
grid systems. The wall normal spacing was varied along the cylinder-flare surface in order to keep the y + 
values within recommended values. This entailed very fine spacing along the face of the flare, 8 micron 
spacing in the region of the cylinder-flare intersection. 

The velocity and thermal boundary layers were measured at ax = —0.06m (just upstream of the cylinder- 
cone corner). Transition was determined to be within the region — 1m < x tr < —0.6m, which includes a 
sizeable portion of the conical nose. Transition was set 0.8m upstream of the flare-cone intersection, and 
this choice did give the best match for the measured velocity and temperature boundary layers. 
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(a) Velocity Boundary Layer 


(b) Enthalpy Boundary Layer 



(c) p/Poo , 20° 


(d) Heat Transfer, 20° 


Figure 6. Kussoy Cylinder Flare Upstream Boundary Layers, Pressure and Heat Transfer, Moo 
Case 
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(a) p/poo, 30° 


(b) Heat Transfer, 30° 




(c) p/poo, 32.5° 


(d) Heat Transfer, 32.5° 



(e) p/poo, 35° ' (f) Heat Transfer, 35° 


Figure 7. Kussoy Cylinder Flare Pressure and Heat Transfer, Moo — 7.05, Separated Cases 
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The 20° case is attached flow. All the models predict pressure well for this case. The heat transfer 
predictions are less stellar, with ail two equation models exhibiting a heat transfer peak not evident in the 
experimental results. The one equation models do not exhibit this overshoot, but underpredict the heating 
rate as the cone surface is downstream of the heating maximum where the other models better predict the 
heat transfer. 

The 30° case is separated flow, The Lag model predicts the separation slightly better that SA, k-cu 
models, but all three are reasonable predictions. The SST model predicts a much larger separation than 
seen experimentally. The heat transfer predictions are in line with the 20° case, with no model performing 
flawlessly, although the SA model’s predictions Eire argueably the best. 

The 32.5° case has more separated flow than the 30° case., The SA, k-cu, and Lag models all predict 
pressure well for generally well for this case, with the Lag model again providing a slightly improved sepa- 
ration prediction. The SST models again gives an overprediction of the extent of the separation. The heat 
transfer predictions in again line with the 20° case, although the overshoot in heating is more pronounced. 

The 35° case not predicted well by any of the models, The models all underpredict the extent of the 
separation, but SST gives the separation prediction closest to the experimental values. The heating rate 
predictions are even less satisfactory, as could be expected given the relatively poor separation predictions 
of the models. 

Various compressibility corrections were attempted, but no universal corrections (in that they reproduced 
or improved earlier predictions, as well as improving these predictions) were found. This remains an ongoing 
area of research for this model. 

C. Overexpanded Nozzle Flow 

This flowfield is an important case for design and analysis of rocket motors utilized over a large range of 
external pressures, such as is the case for space launch vehicles. A ’cold fiow’(non-reacting') rocket nozzle 
expands into an environment 1 with a backpressure sufficiently high to induce separation on the nozzle walls. 
This flowfield is simulated with the help of the chimera capabilities of OVERFLOW. The grid is composed 
of 5 zones: Nozzle Interior, Nozzle Lip, Nozzle Exterior, Near Field Downstream, and Farfield. The nozzle 
interior grid has dimentions 245 x 145, and comprises roughly 3/4 of all the points in the grid system. 

The inflow boundary of the rocket nozzle is accomplished by combining two boundary conditions se- 
quentially, a read of the full conservative conditions followed by a ’nozzle inflow’ condition, which holds p t 
and T t constant, and zeros the radial and circumferential velocities. The static pressure and axial velocity 
are allowed to adjust. The walls are treated as adiabatic no slip surfaces, and the nozzle is immersed in a 
Freestream with M.^ = 0.2, with characteristic boundary conditions. 

Three cases are discussed in this paper, all involving separated flow in the nozzle interior. These cases 
are parameterized by the ratio of the rocket nozzle chamber pressure to the ambient farfield pressure (NPR). 

The least separated flow discussed in this paper has a NPR of 66(Fig. 9(a)). For this case, all of the 
models predict the separation location adequately, and also do a reasonable job of prediction of the pressure 
in the separation region. The post separation presssure is reasonably well predicted. 

A more seperated case has a NPR of 44(Fig. 9(b)), and has more extensive separation. For this case, 
the Lag model did not want to converge on a single location for the separation, and hence was then run in 
time accurate mode using dual timestepping. The time dependent process converged on a roughly constant 
separation/shock location, The range of the separation location for the relaxation includes the position 
predicted by the time accurate simulation. The time accurate solution did still exhibit some unsteadiness, 
and hence the ’steady state’ solution for this case would be suspect in any event. 

The SA, k-cu and SST models converged to a steady state solution without recourse to the time accurate 
simulation, but miss the separation location, and the pressure post separation. The k-cu model consistently 
underpredicts the separation extent, SST overpredicts the separation extent, and all three miss the post 
separation pressure level, which is more accurately predicted by the Lag model. 

The final case has a NPR of 34(Fig. 9(c)), and the most extensive separation. Again, the Lag model 
did not want to converge on a single location for the separation, and hence was then run in time accurate 
mode using dual timestepping. The time dependent process again converged on a separation/shock location, 
though the variation of the shock location and post shock pressure were larger than exhibited by the previous 
case. 
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(b) NPR=44 



(c) NPR=34 

Figure 8. Near Field Mach Number Contours, Lag Model 
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(a) NPR = 66 


(b) NPR = 44 


(c) NPR = 34 


Figure 9. Nozzle Interior Wall Pressure Predictions for Various Backpressures. 


D. Shuttle 



(a) Overset Grid System 


(b) Heating 


Figure 10. Space Shuttle Orbiter, Moo =6.02, Rei_ = 600 x 1 0 3 , <x = 40° 


After examining several axi-symmetric test cases, the lag model was applied to the Space Shuttle Orbiter 
at Moo = 6.02, a = 40°. This configuration represents a complex geometric flowfield at a high angle of 
attack. Freestream conditions were = 62.31°K, p,*, = 2032. 2Pa, = 0.1136kg/m 3 . The model length 

is ,25m Wall conditions were modelled as a constant temperature (T w = 300°K), no slip smooth wall. As 
this is attached flow, the most relevant comparison is the surface heating rate. 

The shuttle data was obtained by Berry and Hamilton 11 in Langley’s 20-inch Mach 6 wind tunnel. 12 It 
is a conventional blowdown facility with a 0.508 x 0.5207m rectangular test section. Electric heaters can 
vary the stagnation temperature up to 590K, and the normal operating stagnation pressure is approximately 
7xl0 5 to 3.7x10 s Pa. Data was obtained by using phosphor thermography techniques. This method uses 
ceramic models that are coated with phosphors that when illuminated with ultraviolet light, fluoresce in 
two regions of the visible spectrum. The fluorescence intensity is dependent on the surface temperature. By 
taking fluorescence intensity images with a color video camera and calibrating the temperature prior to the 
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test, heat transfer can be calculated based on the surface temperature time histories. The thickness of the 
phosphor paint coating is approximately 0.001 inches thick. 

The wind tunnel data was obtained in order to assess the effects of discrete roughness elements on 
transition. The trips were located at an X/L location of 0.258 This allowed for a majority of the windward 
surface to be turbulent. Computations were ran either fully turbulent or fully laminar. 

Grid generation for the Shuttle Orbiter was based on the recommendations given by Chan et al 13 for 
overset grid generation. The orbiter was modelled using six zones. The grid for this case is pictured in 
Figure 10(a). The pitch plane of the orbiter was used as a symmetry plane with y-symmetry. Special care 
was take to ensure each grid had perpendicular grid lines to the pitch plane. Each grid extends from the 
surface to the outer boundary (upstream'of the shock). A minimum of 5 grid points overlap between regions. 

The outer boundary of each grid was fitted to the shock shape for the test condition. Stretching ratios 
are below 1.2 on the surface grids and in volume grids in the off-body direction. Volume grid generation 
took approximately one working day. The total number of grid points for the Shuttle Orbiter grid system 
was nearly 8 million with each grid having 85 points in the off-body direction. Wall spacing values were 
evaluated based on a Re ce ii value given by Gnoffo. 1,1 Re ce ii = pw °" Az < 1 was maintained. The resulting 
y + values were on the order of 0.1 with the chosen wall spacings. 

Heat transfer rates are compared in Figure 10(b). The laminar results compare within the experimental 
uncertainty until the flow was tripped, although they are consistently high. After approximately X/L = 0.258, 
the tripped flow became fully turbulent. All three turbulent predictions are higher than the experimental 
data, although the 2 equation model predictions are consistently higher than the SA model. 

IV. Conclusions 

The Lag model is a suitable model for computing high speed flows, including flows with separation. Skin 
friction is predicted predicted well for attached fiowfields. Good separation predictions are obtained for both 
external and internal geometries. The prediction of the nozzle flow cases using a time accurate calculation 
is exciting. Heat transfer predictions need improvement for separated flows especially. The model has the 
advantages of not requiring wall distance, making it suitable for more advanced computational methods, and 
has shown some promise in computing separated fiowfields when run in a time accurate mode. 
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